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摘要 : 本 文采 用 平衡 分 子 动力 学 方法 研究 了 纳米 流体 中 国体 颗粒 周围 的 吸附 层 厚度 及 结构 ， 并 通过 Green-Kubo 原理 计算 了 吸 
附 层 内 导热 系数 的 分 布 。 研 究 表明 ， 纳 米 颗粒 周围 存在 0.5nm 厚 的 吸附 层 ， 层 内 原子 数 密度 提高 了 50-60%， 导 热 系数 提高 
200-300%。 分 析 发 现 ， 吸 附 层 导热 系数 的 提高 是 由 于 原子 数 密度 及 吸附 层 内 原子 排列 有 序 程度 的 提高 引起 的 。 该 研究 结果 有 望 
为 建立 合理 的 预测 纳米 流体 的 导热 系数 模型 提供 关键 参数 及 理论 指导 。 
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Molecular dynamics study of contribution of the particle surface absorption layer to the 


thermal conductivity of nanofluids 
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Abstract: In this article we studied the thickness and structure of absorption layer around nanoparticle and derived 
the thermal conductivity distribution based on the principle of Green-Kubo by equilibrium molecular dynamics 
simulation. The results showed that there is a 0.5-nm-thickness absorption layer with 50-60% enhancement on atom 
number density and 200-300% enhancement on thermal conductivity. It was found after analysis that the 
enhancement of thermal conductivity(TCE) is the result of both the enhancement of number density and the more 
ordered alignment of atoms within the absorption layer. Our work is expected to provide critical parameters for the 
establishment of the mathematical model for prediction of the thermal conductivity of various nanofluids. 
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0 引言 


纳米 流体 是 以 一 定 的 比例 和 方式 在 液体 中 添加 纳米 级 尺 
度 的 金属 或 非 金 属 固体 颗粒 形成 的 一 类 新 的 传 热 冷却 介质 
巾 。 大 量 的 实验 研究 表明 ， 纳 米 流体 的 导热 系数 由 于 纳米 颗 
粒 的 加 入 而 实现 显著 的 提高 中。 然而 ， 传 统 的 基于 有 效 介质 
理论 的 预测 固 液 两 相 混 合 物 导热 系数 的 经 典 模型 "在 被 用 
于 预测 纳米 流体 的 导热 系数 时 ， 其 预测 结果 往往 较 实 验 值 偏 
小 。 为 了 有 效 预 测 纳米 流体 的 导热 系数 ， 国 内 外 研究 者 基于 
纳米 颗粒 自身 的 特性 及 影响 ， 提 出 了 不 同 的 热传导 机 理 ， 归 
纳 起 来 主要 有 以 下 三 种 : CIO 颗粒 的 布衣 运动 及 其 引起 的 微 
PAEA, (2) 颗粒 的 团聚 (“站 1，(3) 颗粒 表面 的 液体 吸附 
层 [20-25] | 

2001 4E, Choi"? HH aK URLS Se TE ERIN, 1% 
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热 系 数 也 大 于 纯 液 体 。 在 此 理论 
基础 上 ，Yu 和 Choi? 7 e A fi tH T A EJ i Maxwell and 
Hamilton-Crosser #i#!. Leong’), Xie? Al Rizvi "EE A WW 
完工 作 也 充分 考虑 了 颗粒 表面 的 吸附 层 发 展 出 了 不 同 的 预测 
模型 ， 这 些 模 型 与 一 部 分 实验 结果 吻合 的 很 好 ， 但 都 不 具有 
普 适 性 。 
吸附 层 的 厚度 及 导热 系数 是 上 述 模型 的 两 个 重要 参数 ， 
但 这 些 模型 本 身 并 不 能 得 到 这 两 个 参数 。 因 此 ， 有 必要 通过 
其 它 手段 对 吸附 层 进行 研究 。2009 E, Gerardi OR Fl zd 
振 实 验方 法 测量 出 在 氧化 铝 纳米 颗粒 表面 形成 了 厚度 为 1.4 
nm 的 液体 吸附 层 。 然 而， 就 目前 的 报道 而 言 ， 这 样 的 实验 结 
果 是 非常 有 限 的 PC。 而 另 一 方面 ， 在 纳米 级 尺度 下 ， 分 子 动 
力学 的 快速 发 展 为 吸附 层 的 研究 提供 了 一 种 有 效 手 段 睛 汪 。 
本 文采 用 平衡 分 子 动力 学 首次 计算 出 吸附 层 内 导热 系数 的 分 
布 ， 并 对 层 内 原子 数 密度 的 分 布 曲线 及 影响 层 内 的 导热 系数 
的 因素 进行 分 析 研 究 。 

1 模拟 方法 

分 子 动力 学 是 从 微观 层面 研究 纳米 流体 导热 机 理 的 有 效 
手段 CE。 分 子 动力 学 模拟 根据 牛顿 运动 力学 从 每 个 原子 的 初 
始 位 置 及 初始 速度 出 发 ， 由 原子 间 的 势能 函数 求解 出 每 个 原 
子 所 受 的 力 ， 从 而 获得 原子 下 一 时 刻 的 位 置 及 速度 ， 依 次 可 
得 到 系统 内 每 个 原子 随时 间 的 运动 过 程 。 原 子 的 微观 运动 可 
以 用 来 解释 物质 在 宏观 层面 可 观测 的 物理 特性 。 

本 文选 择 氢 (Ar) 原子 构成 基 液 ， 铜 原子 (Cu) 或 银 原 
子 (Ag) 构成 固体 纳米 颗粒 ， 由 于 Ar 原子 的 势能 函数 及 参 
数 与 实验 值 吻合 的 很 好 ， 上 述 纳米 流体 体系 被 各 种 实验 及 理 
论 研究 广泛 使 用 请 23。 


图 1 原子 的 初始 位 置 图 


Fig. 1 The original position of atoms 

如 图 1 所 示 ， 本 文 首先 建立 一 个 边 长 为 10 nm 的 立方 体 
作为 模拟 盒子 ， 在 盒子 的 中 心 划分 出 一 个 球形 区 域 放置 纳米 
颗粒 ， 并 将 纳米 颗粒 限制 在 比 其 半径 大 0.1 nm 的 球形 区 域 以 
方便 后 续 的 导热 系数 计算 。Ar，Cu 和 Ag 原子 均 按照 面 心 立 
方 堆积 CFCC) Jp X HE ps sh dA AY MAE 
573.06,361.49,407.30 pm， 其 中 Ar ECT HE sc i GPL I a 
密度 等 于 1410 kg/m? 计算 而 得 。 

Ar-Ar, Cu-Cu 和 Ar-Cu 原子 间 均 采用 L-J 势能 函数 ， 其 
形式 如 下 : 


U(x) - 0? zm QD 


其 中 , rj 为 原子 i 和 j 之 间 的 距离 ,x 为 截断 半径 ， 在 此 
取 1.5 nm, e fll c 为 势能 函数 的 能 量 参 数 和 长 度 参数 ， 不 同 
类 原子 间 的 势能 参数 可 以 由 Lorentz-Berthelot 混合 法 则 取得 ， 
即 
Sss t On 


Os 2 


£j = Vessel (2) 


本 文 模拟 所 用 到 原子 间 的 势能 参数 由 文献 所 得 ， 具 休 如 下 表 
Bist. 


表 1 Cu/Ar 相互 作用 的 L-J 势能 参数 
Table 1 L-J parameters used for Cu-Ar nanofluids 


原子 对  o/nm 1075] 
Ar-Ar 0.3405 1.654 
Cu-Cu 0.2338 65.581 
Ag-Ag 0.2644 5524 


Cu-AT 
Ag-Ar 


0.2872 10.415 
0.3025 9.559 


模拟 中 采用 周期 性 边界 条 件 ， 时 间 步 长 选取 0.1 fs, KA 
温度 设 定 为 85 k, 体系 在 NVT 系 综 下 运行 2000000 步 ， 再 在 
NVE 系 综 下 运行 2000000 步 ， 经 过 充分 弛 瑰 使 体系 达到 平衡 
状态 ， 记 录 此 时 各 区 域 的 原子 数 ， 最 后 在 NVE 系 综 下 运行 
1000000 步 计 算 导 热 系 数 。 根 据 Green-Kubo 原理 ， 导 热 系数 
可 通过 对 微观 热流 的 自 相关 函数 积分 求 得 ， 即 


1 了 1 7 
K- J.(0)J.(t)dt?- J(0)J(t)dt 
A A frat 
( 3 ) 
1 
Xen 2 (4) 
i i<j 


其 中 ，Ks EBOOK SHR V 是 计算 区 域 的 体积 ，7 是 
体系 的 温度 ,J 是 计算 区 域 原子 的 微观 热流 ，e; 是 原子 i 的 总 
能 量 (势能 与 动能 的 总 和 )，vi 是 原子 i 的 速度 矢量 ,ff 是 原 
子 i 与 j 之 间 的 作用 力 ，xy 是 原子 i 与 j 之 间 的 距离 。 


2 结果 与 讨论 


2.1 平衡 状态 的 确定 
在 计算 导热 系数 之 前 ， 保 证 体系 达到 平衡 状态 是 至 关 习 
要 的 。 当 体系 总 原子 数 较 多 时 (9000 个 原子 )， 不论 是 否 
于 平衡 状态 , 体系 的 温度 或 能 量 都 会 很 快 稳定 在 某 一 定 值 5。 
本 文 的 模拟 体系 中 原子 总 数 超过 20000 个 ， 因 此 ， 仅 仅 通过 
温度 或 能 量 的 涨 落 幅度 已 不 足以 判断 体系 是 否 进入 平衡 状 
态 。 在 有 关 纳 米 流体 吸附 层 的 分 子 动力 学 研究 中 ， 固 体 颗粒 
周围 原子 数 密度 随时 间 的 变化 是 一 种 行 之 有 效 的 平衡 判断 标 
YEP), A 2 可 看 出 ， 在 运行 到 500000 25 (50ps) 时 ， 体 系 
的 温度 波动 范围 已 经 很 小 ， 我 们 选择 从 此 刻 开始 进行 分 析 。 
如 图 3 所 示 ， 横 坐标 表示 球 壳 区 域 与 固体 纳米 颗粒 表面 
的 径 向 距离 ， 纵 坐标 表示 该 区 域 的 原子 数 密度 n, BREK 
域内 的 原子 数 Y 及 体积 V Sp ET A: 
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n=— (5) 
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步 数 图 2 ath 
豫 过 程 中 体系 的 温度 变化 


Fig. 2 The evolution oftemperature during relaxation 
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不 同时 刻 液体 原子 数 密度 比较 
Fig. 3 Comparison of number density at different simulation 
times 
由 图 3 可 以 看 出 ，50ps 与 100ps 两 时 刻 的 原子 数 密度 分 
布 相 比较 ， 靠 近 固 体 纳 米 颗 粒 表面 的 区 域 原子 数 密度 大 幅度 
上 升 ， 说 明 在 这 段 时 间 内 体系 中 的 原子 由 于 固 液 原子 间 的 相 
互 作用 向 固体 纳米 颗粒 表面 迁移 ， 体 系 并 未 进入 平衡 状态 。 
100ps 与 200，300，400ps 相 比 ， 原 子 数 密度 分 布 波动 很 小 ， 
可 以 判断 出 400ps 时 体系 已 经 达到 平衡 状态 。 
通过 弛 驳 过 程 使 体系 进入 平衡 状态 , 就 可 以 在 NVE 系 综 
下 进行 导热 系数 的 计算 ， 而 公式 (3) 中 的 积分 时 间 原则 上 取 无 
限 长 ， 在 实际 计算 过 程 中 计算 时 间 不 可 能 无 限 长 ， 但 积分 时 
间 选 取 过 短 ， 计 算 误 


图 4 导热 
- herd 随 积分 时 
E " a >, 
Ne 省 ， 间 的 变化 
Up " Fig. 4 
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thermal conductivity and integral time 
差 会 偏 大 。 采 用 平衡 分 子 动力 学 的 方法 计算 导热 系数 ， 可 以 
根据 热流 密度 或 导热 系数 的 计算 结果 随 积 分 时 间 的 变化 选取 
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合理 的 积分 时 间 PU. WE 4 所 示 ， 可 以 看 出 随 着 积分 时 间 变 
长 ， 导 热 系数 逐渐 趋 近 于 某 一 常数 ， 计 算 误差 趋 近 于 0。 本 
模拟 选取 100000 27 (100ps) 作为 积分 时 间 是 合理 的 。 
2.2 原子 数 密度 的 分 布 

在 分 子 动力 学 模拟 中 ， 纳 米 颗 粒 周 围 原 子 数 密度 的 提高 
被 用 来 研究 吸附 层 "*”。 本 文 将 纳米 颗粒 周围 的 区 域 分 成 一 
系列 紧邻 的 厚度 为 0.1 nm 的 球 觉 区域， 计算 每 个 区 域 的 原子 
数 密度 ， 得 到 原子 数 密度 分 布 图 。 如 图 5， 可 以 看 出 在 
0.2,0.5,0.8 nm 位 置 ， 原 子 数 密度 出 现 峰 值 ， 在 0.8 nm 之 外 原 
子 数 密度 分 布 均匀 可 视 为 纯 液 氢 的 原子 数 密度 。 相 比 于 纯 液 
氢 的 原子 数 密度 21/nm*， 在 0-0.2 nm 区 域 平均 原子 数 密度 为 
46mmm，0-0.5 nm 区 域 平均 原子 数 密度 为 29/mm ， 这 表明 纳 
米 颗 粒 周围 形成 了 厚度 约 为 0.5 nm 的 吸附 层 。 平 均 而 言 ， 厚 
度 为 0.5 nm 的 吸附 层 原 子 数 密度 提高 了 50-60%. 

同时 ， 从 图 5 中 也 可 以 看 出 ,在 0.1，0.3 nm 等 处 层 内 原 
子 数 密度 明显 小 于 纯 液 毛 的 原子 数 密度 ， 表 明 吸 附 层 内 原子 
数 密 度 的 分 布 是 极 不 均匀 的 ， 其 原因 是 吸附 层 结构 的 形成 不 
仅 取决 于 固 液 间 相互 作用 ， 液 体 间 的 相互 作用 也 是 不 可 忽视 
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径 向 距离 am) 
图 7 银 纳米 颗粒 周围 原子 数 密度 的 分 布 曲线 


Fig.7 Number density distribution curve around Ag nanoparticle 

原子 数 密度 ， 但 实际 曲线 表现 为 先 急剧 的 减 小 再 增加 到 
另 一 个 峰值 ， 这 是 因为 在 0.2-0.5 nm 区 域内 ， 液 体 间 的 相互 
作用 占据 主导 地 位 ， 由 于 0.1-0.2 nm 区 域内 原子 数 密度 极 大 ， 
这 对 相 邻 区 域内 原子 的 作用 力也 得 到 了 加 强 ， 对 距离 较 近 的 
0.2-0.3 nm 区 域内 的 原子 表现 出 斥 力作 用 使 其 原子 数 密度 急 
剧 减 小 ， 对 距离 较 远 的 0.4-0.5 nm 区 域内 的 原子 引力 增加 使 


的 。 结 合 图 6， 在 0-0.1 nm 区 域 ， 由 于 距离 太 近 ， 固 液 间 表 
现 为 斥 力作 用 ， 因 此 原子 数 密度 较 小 ;在 0.1-0.2 nm 区 域 ， 

由 于 固 液 间 的 相互 吸引 作用 ， 原 子 数 密度 出 现 极 大 值 ;在 
0.2-0.5 nm 区 域 ， 如 果 只 考虑 固 液 间 的 相互 作用 ,原子 数 密度 
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图 5 铜 纳米 颗粒 周围 原子 数 密度 的 分 布 曲 线 


Fig. 5 Number density distribution curve around Cu nanoparticle 
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图 6 两 原子 间 的 势能 随 距离 的 变化 


Fig. 6 Evolution ofpotential energy with distance between two 


atoms 


其 原子 数 密度 增 大 。0.5 nm 之 后 固 液 间 相互 作用 可 以 忽略 
原子 数 密度 出 现 的 极 大 极 小 值 都 可 以 依次 解释 。 王 楠 等 人 5 
的 研究 工作 表明 吸附 层 呈 现 出 双 层 结构 ， 类 似 于 扩散 双 电 层 
模型 ， 内 层 厚 度 约 为 0.42nm， 内 外 层 总 厚度 约 为 0.76nm， 这 
与 我 们 的 研究 结果 吻合 ， 但 对 出 现 该 结构 的 深层 次 原因 ， 该 
文 并 未 进一步 深入 探讨 。Cui 等 人 "的 研究 中 ， 随 时 间 变 化 ， 
原子 数 密度 曲线 整体 下 降 ， 表 明 模 拟 过 程 中 体系 的 总 原子 数 
在 减少 ， 就 其 采用 的 周期 性 边界 条 件 而 言 ， 出 现 这 样 的 结果 
是 不 合理 的 。 

随 着 纳米 颗粒 直径 从 1 nm 增 大 到 3 nm， 吸 附 层 的 厚度 
及 结构 并 未 出 现 明显 变化 。 对 比 图 7， 可 看 出 不 同 材料 的 纳 
米 颗 粒 ， 吸 附 层 的 厚度 及 结构 变化 不 大 。 
2.3 导热 系数 的 分 布 

本 研究 采用 平衡 分 子 动力 学 方法 , 基于 Green-Kubo 原理 ， 
在 体系 达到 平衡 状态 后 ， 计 算 每 个 球 壳 区 域 的 导热 系数 ， 可 
以 得 到 导热 系数 的 分 布 曲线 (不 包含 0-0.1 nm 区 域 的 导热 系 
数 )， 结 果 如 图 8 RAO 所 示 。 平均 而 言 ， 在 0.1-0.2 nm 区 域 
内 ， 导 热 系 数 提高 了 300%-600%; 在 0.1-0.5 区 域内 ， 导 热 系 
数 提高 了 约 200%-300%。 这 表明 吸附 层 内 的 导热 系数 大 于 纯 
RiT o 
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力学 方法 计算 纳米 流体 的 导热 系数 ， 两 种 不 同 的 方法 结果 是 
一 致 的 ， 相 差 值 是 微小 到 可 忽视 的 。 而 Liang APRA 4E 
平衡 分 子 动力 学 方法 ， 对 于 Cu/Ar 体系 的 研究 结果 表明 ， 吸 
附 层 的 导热 系数 提高 了 约 60%， 这 与 我 们 的 结果 差距 较 大 ， 
其 原因 主要 是 其 吸附 层 的 厚度 选取 为 1.158nm。 当 我 们 的 吸 
附 层 厚度 取 1.1nm 时 ， 吸 附 层 导热 系数 提高 了 40%-50%， 则 


小 于 60%。 这 主要 是 由 于 Liang 的 模拟 体系 中 ， 吸 附 层 形成 
于 铜 原 子 组 成 的 墙壁 而 不 是 球形 颗粒 ， 吸 附 层 结构 是 矩形 层 
状 的 而 不 是 球 壳 状 。 这 就 意味 着 相对 于 我 们 的 模型 ，Liang 等 
人 的 模型 中 具有 高 导热 系数 的 0.1-0.2 nm 薄 层 的 体积 有 显著 
提高 ， 从 而 解释 了 其 结果 与 我 们 结果 的 差 值 。 

另外 ,对 比 原子 数 密度 与 导热 系数 分 布 图 , 可 看 出 , 图 $ 
和 图 7 中 的 曲线 分 别 与 图 8 和 图 9 中 对 应 的 曲线 变化 趋势 同 
步 一 致 。 在 0.2，0.5，0.8 nm (甚至 包括 1.1 nm) 处 原子 数 密 
度 出 现 峰 值 , 导热 系数 也 出 现 相 应 的 峰值 ; 在 0.3, 0.6, 0.9 nm 
处 原子 数 密度 出 现 极 小 值 ， 导 热 系数 也 出 现 极 小 值 ， 同 时 ， 
也 可 以 看 出 ， 任 意 两 区 域 相 比 ， 原 子 数 密度 大 ， 则 导热 系数 
也 大 ， 反 之 亦 然 。 这 表明 导热 系数 的 提高 与 原子 数 密度 密切 
相关 。 


Ed —e— ^nm 
€ —e— 2nm 
E 1 —e— 3nm 
E 0.44 
ES e 
K 0.24 
OP :8 6:9 * 0:9. 9-09.9.9:9 9.90 
径 向 距离 (nm) 
图 8 铜 纳米 颗粒 周围 导热 系数 的 分 布 曲线 


Fig. 8 thermal conductivity distribution curve around Cu 
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图 9 银 纳 米 颗粒 周围 导热 系数 的 分 布 曲线 


Fig. 9 thermal conductivity distribution curve around Ag 
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纳米 颗粒 周围 单个 Ar 原子 导热 系数 的 分 布 曲线 


Fig. 10 thermal conductivity distribution curve of serial Ar atom 


around Cu nanoparticle 
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图 11 银 纳米 颗粒 周围 单个 Ar 原子 导热 系数 的 分 布 曲线 


Fig. 11 thermal conductivity distribution curve of serial Ar atom 


around Ag nanoparticle 
为 了 探究 原子 数 密度 与 导热 系数 的 关系 ， 我 们 定义 了 如 
下 的 无 量 纲 原子 数 密度 n: 


n-- (6) 


Eh, n 为 各 球 壳 区 域 的 原子 数 密度 ，mo 为 纯 液 氮 的 原子 
数 密度 ， 取 21l/nm 。 本 文采 用 导热 系数 与 无 量 纲 原子 数 密度 
的 比值 表征 每 个 原子 对 导热 系数 的 影响 ， 可 得 到 图 10 和 图 
11。 可 以 看 出 ， 在 0.1-0.5 nm 的 区 域 ， 每 个 吸附 层 原子 对 导 
热 系 数 的 贡献 也 有 显著 提高 ， 平 均 而 言 提高 了 约 40-50%。 

通过 以 上 分 析 可 看 出 ， 吸 附 层 的 形成 及 其 高 导热 特性 根 
本 原因 在 于 较 强 的 固 液 间 的 相互 作用 。 由 于 固 液 间 的 这 种 相 
互 作用 ， 一 方面 使 其 周围 原子 数 密度 增 大 ;， 另 一 方面 原子 排 
列 的 有 序 化 程度 提高 ， 导 热 结构 得 到 优化 ， 表 现 为 每 个 原子 
对 于 导热 系数 的 贡献 也 有 所 提高 。 


本 文通 过 平衡 分 子 动力 学 方法 详细 研究 了 吸附 层 原子 数 


密度 及 导热 系数 的 分 布 。 主 要 得 到 以 下 结论 : 

C1) 本 文 模拟 表明 吸附 层 的 厚度 约 为 0.5 nm。 相 比 于 纯 
液体 ， 吸 附 层 内 的 原子 数 密度 提高 了 50-60%， 导 热 系 数 提高 
了 200-300%。 该 结果 为 准确 预测 纳米 流体 的 导热 系数 提供 了 
Hed. 

(2) 由 于 固 液 和 液 液 间 原子 的 相互 作用 ， 吸 附 层 内 原本 
数 密度 及 导热 系数 分 布 极 不 均匀 ， 且 在 本 文 的 模拟 中 ， 吸 附 
层 的 厚度 并 未 随 着 纳米 颗粒 粒 径 的 变化 而 出 现 明显 变化 。 

(3) 吸附 层 导热 系数 的 提高 一 方面 是 由 于 原子 数 密度 的 
提高 ， 另 一 方面 是 由 于 原子 排列 有 序 化 程度 提高 ， 两 者 的 影 
响 基本 相当 。 

(4) 吸附 层 形成 是 由 于 固 液 间 的 相互 作用 力 强 于 液体 内 
部 作用 力 ， 这 启示 我 们 在 实际 选取 纳米 流体 时 ， 固 液 间 的 相 
互 作用 也 是 一 个 值得 考虑 的 因素 。 
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